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Role of nonlinear anisotropic damping in the magnetization dynamics of topological solitons 
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The consequences of nonlinear anisotropic damping, driven by the presence of Rashba spin-orbit coupling 
in thin ferromagnetic metals, are examined for the dynamics of topological magnetic solitons such as domain 
walls, vortices, and skyrmions. The damping is found to affect Bloch and Neel walls differently in the steady 
state regime below Walker breakdown and leads to a monotonic increase in the wall velocity above this transition 
for large values of the Rashba coefficient. For vortices and skyrmions, a generalization of the damping tensor 
within the Thiele formalism is presented. It is found that chiral components of the damping affect vortex- and 
hedgehog-like skyrmions in different ways, but the dominant effect is an overall increase in the viscous-like 
damping. 


PACS numbers: 75.60.Ch, 75.70.Kw, 75.75.-c, 75.78.Fg 

I. INTRODUCTION 

Dissipation in magnetization dynamics is a longstanding 
problem in magnetism [1-3]. For strong ferromagnets such as 
cobalt, iron, nickel, and their alloys, a widely used theoretical 
approach to describe damping involves a local viscous form 
due to Gilbert for the Landau-Lifshitz equation of motion, 

dm dm 

— = -yomxHeff + a 0 mx —, (1) 

at at 

which appears as the second term on the right hand side, pro¬ 
portional to the damping constant ao- This equation describes 
the damped magnetization precession about a local effective 
field H e ff = -(1 /i-inM s )6U /dm, which is given by a variational 
derivative of the magnetic energy U with respect to the mag¬ 
netization field described by the unit vector m, with yo = po7 
being the gyromagnetic constant and M s is the saturation mag¬ 
netization. Despite the multitude of physical processes that 
underlie dissipation in such materials, such as the scattering 
of magnons with electrons, phonons, and other magnons, the 
form in Eq. (1) has proven to be remarkably useful for describ¬ 
ing a wide range of dynamical phenomena from ferromagnetic 
resonance to domain wall motion. 

One feature of the dissipative dynamics described in Eq. (1) 
is that it is local, i.e., the damping torque only depends on the 
local magnetization and its time dependence. With the ad¬ 
vent of magnetic heterostructures, however, this restriction of 
locality has been shown to be inadequate for systems such 
as metallic multilayers in which nonlocal processes can be 
important [4], A striking example involves spin pumping, 
which describes how spin angular momentum can be dissi¬ 
pated in adjacent magnetic or normal metal layers through the 
absorption of spin currents generated by a precessing magne¬ 
tization [5, 6]. Early experimental observations of this phe¬ 
nomena involved iron films sandwiched by silver layers [7] 
and permalloy films in close proximity with strong spin-orbit 
normal metals such as palladium and platinum [8, 9], where 
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ferromagnetic resonance line widths were shown to depend 
strong on the composition and thickness of the adjacent lay¬ 
ers. Such observations also spurred other studies involving 
ferromagnetic multilayers separated by normal metal spacers, 
where spin pumping effects can lead to a dynamic coupling 
between the magnetization in different layers [10, 11]. In 
the context of damping, such dynamic coupling was shown to 
give rise to a configuration dependent damping in spin-valve 
structures [12, 13]. 

A generalization of the spin-pumping picture in the context 
of dissipation was given by Zhang and Zhang, who proposed 
that spin currents generated within the ferromagnetic material 
itself can lead to an additional contribution to the damping, 
provided that large magnetization gradients are present [14]. 
This theory is based on an sd model in which the local mo¬ 
ments ( Ad) are exchange coupled to the delocalized conduc¬ 
tion electrons (3s), which are treated as a free electron gas. 
The spin current “pumped” at one point in the material by 
the precessing local moments are dissipated at another if the 
current encounters strong spatial variations in the magneti¬ 
zation such as domain walls or vortices - a mechanism that 
can be thought of as the reciprocal process of current-induced 
spin torques in magnetic textures [15-18], For this reason, 
the mechanism is referred to as “feedback” damping since the 
pumped spin currents generated feed back into the magnetiza¬ 
tion dynamics in the form of a dissipative torque. This addi¬ 
tional contribution is predicted to be both nonlinear and non¬ 
local, and can have profound consequences for the dynamics 
of topological solitons such as domain walls and vortices as a 
result of the spatial gradients involved. Indeed, recent experi¬ 
ments on vortex wall motion in permalloy stripes indicate that 
such nonlinear contributions can be significant and be of the 
same order of magnitude as the usual Gilbert damping char¬ 
acterized by oro in Eq. (1) [19]. 

An extension to this feedback damping idea was proposed 
recently by Kim and coworkers, who considered spin pump¬ 
ing involving a conduction electron system with a Rashba 
spin-orbit coupling (RSOC) [20]. By building upon the 
Zhang-Zhang formalism, it was shown that the feedback 
damping can be expressed as a generalization of the Landau- 
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Lifshitz equation [14, 20], 

3m <9m 

— = -y 0 m x H eff + m x 2) LL (m) • —, (2) 

at at 

where the 3x3 matrix £>ll represents the generalized damping 
tensor, which can be expressed as [ 20 ] 

®ll = a oS u + tj ^ ( F ki + a R e 3fa ) (fy + a R e 3 ^. (3) 

k 

Here, a o is the usual Gilbert damping constant, q = 
glu B f>G Q /(4e 1 M s ) is a constant related to the conductivity Go 
of the spin bands [14], = (dm/dx^j are components of 

the spatial magnetization gradient, a R = la^jnJTi 2 is the 
scaled Rashba coefficient, is the Levi-Civita symbol, and 
the indices ( ijk ) represent the components (xyz) in Cartesian 
coordinates. In addition to the nonlinearity present in the 
Zhang-Zhang picture, the inclusion of the a R term results 
in an anisotropic contribution that is related to the underly¬ 
ing symmetry of the Rashba interaction. Numerical estimates 
based on realistic parameters suggest that the Rashba con¬ 
tribution can be much larger than the nonlinear contribution 
q alone [ 20 ], which may have wide implications for soliton 
dynamics in ultrathin ferromagnetic films with perpendicular 
magnetic anisotropy, such as Pt/Co material systems, in which 
large spin-orbit effects are known to be present. 

In this article, we explore theoretically the consequences 
of the nonlinear anisotropic damping given in Eq. (3) on the 
dynamics of topological magnetic solitons, namely domain 
walls, vortices, and skyrmions, in which spatial gradients can 
involve 180° rotation of the magnetization vector over length 
scales of 10 nm. In particular, we examine the role of chiral¬ 
ity in the Rashba-induced contributions to the damping, which 
are found to affect chiral solitons in different ways. This ar¬ 
ticle is organized as follows. In Section II, we discuss the 
effects of nonlinear anisotropic damping on the dynamics of 
Bloch and Neel domain walls, where the latter is stabilized 
by the Dzyaloshinskii-Moriya interaction. In Section III, we 
examine the consequences of this damping for vortices and 
skyrmions, and we derive a generalization to the damping 
dyadic appearing in the Thiele equation of motion. Finally, 
we present some discussion and concluding remarks in Sec¬ 
tion IV. 


II. BLOCH AND NEEL DOMAIN WALLS 

The focus of this section are domain walls in ultrathin 
films with perpendicular magnetic anisotropy. Consider a 
180° domain wall representing a boundary separating two 
oppositely magnetized domains along the x axis, with z be¬ 
ing the uniaxial anisotropy axis that is perpendicular to the 
film plane. We assume that the magnetization remains uni¬ 
form along the y axis. The unit magnetization vector m(x, 1) 
can be parametrized in spherical coordinates ( 9 , 0 ), such that 
m = (sin 0 cos 0, sin 0 sin 0, cosfl). With this definition, the 


spherical angles for the domain wall profile can be written as 

6(x, t) = 2 tan -1 exp |±- ^ j , 

4Kx , t) = 4>o(t), (4) 

where X f) (t) denotes the position of the domain wall, A = 
yjA/Ko represents the wall width parameter that depends on 
the exchange constant A and the effective uniaxial anisotropy 
K{), and the azimuthal angle (Mi) is a dynamic variable but 
spatially uniform. The anisotropy constant, K {) = K„ - 
l-i(_)M 2 J 2, involves the difference between the magnetocrys¬ 
talline ( K „) and shape anisotropies relevant for an ultrathin 
film. In this coordinate system, a static Bloch wall is given by 
<po = ±7:/ 2, while a static Neel wall is given by 0o — 0, n. A 
positive sign in the argument of the exponential function for 
9 in Eq. (4) describes an up-to-down domain wall profile go¬ 
ing along the +x direction, while a negative sign represents a 
down-to-up wall. 

To determine the role of the nonlinear anisotropic damping 
term in Eq. (3) on the wall dynamics, it is convenient to com¬ 
pute the dissipation function W(X o, (bo) for the wall variables, 
where the notation Xq = d,Xo, etc., denotes a time derivative. 
The dissipation function per unit surface area is given by 

m r°° 

w (x 0 , <Po) = zr 2 I dx riii £>^ L (m) ./’ ( 5 ) 

u—co 

where m, = mi [x - Xo(t), (po(t)] and the Einstein summation 
convention is assumed. By using the domain wall ansatz 
(4), the integral in Eq. (5) can be evaluated exactly to give 
VV = VV'o + Wnl. where Wo represents the usual (linear) Gilbert 
damping. 


Wq = a 0 
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while Wnl is the additional contribution from the nonlinear 
anisotropic damping, 


Wnt = 


M, A 


1 X 2 

-a 3 sin 2 0 o (f)^| 


+ I Aqti ± -a 2 cos M t) + (*3 cos 2 0 o( 0 112 


(7) 


where a i = / 7 /A 2 , 0-2 = //Ar/A, and a 3 = // a 2 are dimen¬ 
sionless nonlinear damping constants. In contrast to the linear 
case, the nonlinear anisotropic dissipation function exhibits a 
configuration-dependent dissipation rate where the prefactors 
of the X 2 and 0^ terms depend explicitly on 0o(f). 

In addition to the nonlinearity a chiral damping term, pro¬ 
portional to a 2 , appears as a result of the Rashba interaction 
and is linear in the Rashba coefficient a R . The sign of this 
term depends on the sign chosen for the polar angle 9 in the 
wall profile (4). To illustrate the chiral nature of this term, we 
consider small fluctuations about the static configuration by 
writing 0 o(f) = 0 o + 6</>(t), where (50(f) <s 7 r is a small angle. 
This approximation is useful for the steady state regime below 
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Walker breakdown. For up-to-down Bloch walls (0o = +n/2), 
the nonlinear part of the dissipation function to first order in 
6<p(t ) becomes 


M s A 

Wnl,B loch ~ - 

y 


03 ^0 
3 A 2 




v2 


( 8 ) 


The quantity C, = +1 is a component of the chirality vec¬ 
tor [21], 

1 f°° 

C = - I dx m x <9 v in, (9) 

— oo 


which characterizes the handedness of the domain wall. For 
a right-handed Bloch wall, 0o = —n/2 and the only nonva¬ 
nishing component is C x = 1, while for a left-handed wall 
(00 = -n/2) the corresponding value is C x — -1. Thus, the 
term proportional to a 2 depends explicitly on the wall chiral¬ 
ity. Similarly for up-to-down Neel walls, the same lineariza¬ 
tion about the static wall profile leads to 


M s A / 2a i „ nao 

IkNL.Neel ~ - —; 1" C y — -1- ff 3 

y \ 3 2 


-.2 


( 10 ) 


where C v = 1 for a right-handed Neel wall (0o = 0) and 
C y = -1 for its left-handed counterpart (0o = zr). Since the 
fluctuation 6cj){t) is taken to be small, the chiral damping term 
is more pronounced for Neel walls in the steady-state velocity 
regime since it does not depend on the fluctuation amplitude 
<J0(f) as in the case of Bloch walls. 

To better appreciate the magnitude of the chirality- 
dependent damping term, it is instructive to estimate numer¬ 
ically the relative magnitudes of the nonlinear damping con¬ 
stants ai,a 2 , Qf 3 . Following [Ref. 20], we assume q = 0.2 nm 2 
and ff R = 10 ~ 10 eV m. If we suppose A = 10 nm, which is 
consistent with anisotropy values measured in ultrathin films 
with perpendicular anisotropy [ 22 ], the damping constants can 
be evaluated to be a\ — 0.002, a 2 = 0.052, and a^ = 1.37. 
Since ao varies between 0.01-0.02 [23] and 0.1-0.3 [24] de¬ 
pending on the material system, the chiral term is compa¬ 
rable to Gilbert damping in magnitude, but remains almost an 
order of magnitude smaller than the nonlinear component a 3 
that provides the dominant contribution to the overall damp¬ 
ing. 

The full equations of motion for the domain wall dynam¬ 
ics can be obtained using a Lagrangian formalism that ac¬ 
counts for the dissipation given by W [25, 26]. For the sake 
of simplicity, we will focus on wall motion driven by mag¬ 
netic fields alone, where a spatially-uniform magnetic field 
H. is applied along the +z direction. In addition, we include 
the Dzyaloshinskii-Moriya interaction appropriate for the ge¬ 
ometry considered [27, 28] when considering the dynamics 
of Neel walls. From the Euler-Lagrange equations with the 
Rayleigh dissipation function, 

^-|^ + ^= 0 , (id 

dt fiX{) dX{) ()X() 

with an analogous expression for 00 , the equations of motion 
for the wall coordinates are found to be 

00 + (<*0 + y sin 2 0o j y = y 0 H z , (12) 




- 0.15 ... 

0 0.1 0.2 0.3 0.4 0 0.1 0.2 0.3 0.4 

a R (eV nm) a R (eV nm) 


FIG. 1. (Color online) Bloch wall dynamics, (a) Steady-state domain 
wall velocity, (X 0 ), as a function of perpendicular applied magnetic 
field, /-IqH-, for several values of the Rashba coefficient, a R . The 
horizontal dashed line indicates the Walker velocity and the arrows 
indicate the Walker transition, (b) The ratio between the Walker ve¬ 
locity, vw, to its linear damping value, vw.o, as a function of a R . (c) 
Deviation in the wall angle from rest at the Walker velocity, d</>w, as 
a function of a R 


X 0 , a. 
A “l a ° + 


2cri na 2 

- 1 -- cos < 
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, + <23 COS 0() | 0o 

2 K ± 


2 //0M.5A qiQM s 


cos 0o sin 0o, (13) 


where D ex is the Dzyaloshinskii-Moriya constant [28] and 
K l represents a hard-axis anisotropy that results from vol¬ 
ume dipolar charges. The Dzyaloshinskii-Moriya interaction 
(DMI) is present in ultrathin films in contact with a strong 
spin-orbit coupling material [29, 30] and favors a Neel-type 
wall profile [31, 32]. The DMI itself can appear as a con¬ 
sequence of the Rashba interaction and therefore its inclu¬ 
sion here is consistent with the nonlinear anisotropic damping 
terms used [20, 33, 34], 

Results from numerical integration of these equations of 
motion for Bloch and Neel walls are presented in Figs. 1 and 
2. We used parameters consistent with ultrathin films with 
perpendicular anisotropy, namely a<) = 0.1, M s - 1 MA/m, 
A = 10 nm, and K ± = ji{yN x M\j2 with the demagnetiza¬ 
tion factor N x = 0.02 [28], To study the dynamics of the 
Dzyaloshinskii (Neel) wall we assumed a value of D ex = 1 
mJ/m 2 , which is much stronger than the volume dipolar in¬ 
teraction represented by K ± and is of the same order of mag¬ 
nitude as values determined by Brillouin light spectroscopy 
in Pt/Co/AhOi films [35]. As in the discussion on numeri¬ 
cal estimates above, we assumed q = 0.2 nm 2 but considered 
several different values for the Rashba coefficient o-r. The 
steady-state domain wall velocity, ( Xq ), was computed as a 
function of the perpendicular applied magnetic field, H z . In 
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the precessional regime above Walker breakdown in which 
< f>o(t ) becomes a periodic function in time, (Xq) is computed 
by averaging the wall displacement over few hundred periods 
of precession. 

For the Bloch case [Fig. 1(a)], the Walker field is observed 
to increase with the Rashba coefficient, which is consistent 
with the overall increase in damping experienced by the do¬ 
main wall. However, there are two features that differ qual¬ 
itatively from the behavior with linear damping. First, the 
Walker velocity is not attained for finite a r, where the peak 
velocity at the Walker transition is below the value reached 
for o' R = 0. This is shown in more detail in Fig. 1(b), where 
the ratio between the Walker velocity, vw, and its linear damp¬ 
ing value, vw.o, is shown as a function of <t- R . The Walker limit 
is set by the largest extent to which the wall angle 0 o can de¬ 
viate from its equilibrium value, 0 o jeq . By assuming 0 = 0 
in the linear regime, we can determine this limit by rearrang¬ 
ing Eqs. 12 and 13 to obtain the following relationship for the 
Bloch wall, 

= - (a 0 + y sin 2 0 O ) sin 20o- (14) 

The angle 0o = 0w for which the right hand side of this 
equation is an extremum determines the Walker limit. In 
Fig. 1(c), we present this limit in terms of the deviation an¬ 
gle, <50w = |0w - 0o,eqli which is shown as a function of o-r. 
As the Rashba parameter is increased, the maximum wall tilt 
possible in the linear regime decreases from the linear damp¬ 
ing value of 7 r/ 4 , which results in an overall reduction in the 
Walker velocity. Second, the field dependence of the wall ve¬ 
locity below Walker breakdown is nonlinear and exhibits a 
slight convex curvature, which becomes more pronounced as 
ctr increases. This curvature can be understood by examining 
the wall mobility under fields, which can be deduced from Eq. 
( 12 ) by setting 0 = 0 , 

=- y °f . 2 H z . (15) 

ao + (as/3) sin” 0 o 

Since the angle 0 O for Bloch walls varies from its rest value of 
0o,eq = ±tt/ 2 at zero field to 0w at the Walker field, the sin 2 0o 
term in the denominator decreases from its maximum value of 
sin 2 0 o, eq = 1 at rest with increasing applied field and therefore 
an increase in the mobility is seen as H. increases, resulting 
in the convex shape of the velocity versus field relation below 
Walker breakdown. 

It is interesting to note that the nonlinear damping terms 
affect the Dzyaloshinskii (Neel) wall motion differently. In 
contrast to the Bloch case, the Walker velocity for increasing 
ctr slightly exceeds the linear damping value, which can be 
seen by the arrows marking the Walker transition in Fig. 2(a) 
and in detail in Fig. 2(b). In addition, the field dependence of 
the velocity exhibits a concave curvature below breakdown, 
which can also be understood from Eq. (15) by considering 
that 0 o instead deviates from the rest value of 0 o, eq = 0 or n 
at zero field. As for the Bloch wall case, the deviation angle 
at breakdown is determined by the value of 0 o that gives an 





FIG. 2. (Color online) Dzyaloshinskii (Neel) wall dynamics, (a) 
Steady-state domain wall velocity, (X 0 ), as a function of perpendic¬ 
ular applied magnetic field, for several values of the Rashba 

coefficient, ur. The horizontal dashed line indicates the Walker ve¬ 
locity and the arrows indicate the Walker transition, (b) The ratio 
between the Walker velocity, vw, to its linear damping value, v w ,o, as 
a function of a R . (c) The wall angle at the Walker velocity, 0 W , as a 
function of a R 

extremum for the right hand side of 

2 H~ { as o \( 7rD ex \ 

-— = -a ( ) + — sin” 0 o -— cos 0 o + sin 20 O , (16) 

N X M S V 3 T )\2K ± A r ) 

and is also seen to decrease with increasing Rashba coeffi¬ 
cient [Fig. 2(c)]. In contrast to the Bloch wall case, how¬ 
ever, changes in 0 w have a comparatively modest effect on the 
Walker velocity. The shape of the velocity versus field curve 
is consistent with experimental reports of field-driven domain 
wall motion in the Pt/Co (0.6 nmVAFOs system [36], which 
possess a large DMI value [35] and harbors Neel-type domain 
wall profiles at equilibrium [37]. 

As the preceding discussion shows, the differences in the 
field dependence of the wall velocity for the two profiles are 
a result of the DMI, rather than the chiral damping term that 
is proportional to as- This was verified by setting as = 0 
for the Neel wall case with D + 0, which did not modify 
the overall behavior of the field dependence of the velocity. In 
the one-dimensional approximation for the wall dynamics, the 
DMI enters the equations of motion like an effective magnetic 
field along the x axis, which stabilizes the wall structure by 
minimizing deviations in the wall angle 0 o(f). 

m. VORTICES AND SKYRMIONS 

The focus of this section is on the dissipative dynamics 
of two-dimensional topological solitons such as vortices and 
skyrmions. The equilibrium magnetization profile for these 
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micromagnetic objects are described by a nonlinear differ¬ 
ential equation similar to the sine-Gordon equation, where 
the dispersive exchange interaction is compensated by dipo¬ 
lar interactions for vortices [38, 39] and an additional uniax¬ 
ial anisotropy for skyrmions [40]. The topology of vortices 
and skyrmions can be characterized by the skyrmion winding 
number Q, 


Q = — jj dxdy m • (3 v m x 0 v m). 


(17) 


While the skyrmion number for vortices (Q = ±1/2) and 
skyrmions (Q = ±1) are different, their dynamics are quali¬ 
tatively similar and can be described using the same formal¬ 
ism. For this reason, vortices and skyrmions will be treated 
on equal footing in what follows and distinctions between the 
two will only be drawn on the numerical values of the damp¬ 
ing parameters. 

A key approximation used for describing vortex or 
skyrmion dynamics is the rigid core assumption, where it is 
assumed that the spin structure of the soliton remains unper¬ 
turbed from its equilibrium state during motion. Within this 
approximation, the dynamics is given entirely by the position 
of the core in the film plane, Xo(f) = [Xo(0, fo(0], which al¬ 
lows the unit magnetization vector to be parametrized as 


d(x,y,t) = 0 O [||x - XqCOII] . 


tfi(x, y,t) - q tan 


y ~ Yp(t) 

x-Xo{t)\ +C 2' 


(18) 


where q is a topological charge and c is the chirality. An il¬ 
lustration of the magnetization field given by the azimuthal 
angle cf>(x,y) is presented in Fig. 3. q = 1 corresponds to a 
vortex or skyrmion, while q = -1 represents the antivortex or 
antiskyrmion. 

The dynamics of a vortex or skyrmion in the rigid core ap¬ 
proximation is given by the Thiele equation, 

C\T 7 

GxXo+£>t-X 0 = - — , (19) 

<3X 0 

where 

G = — JJ dxdy sin(0) (V0 x V0) (20) 

is the gyrovector and f/(Xo) is the effective potential that is ob¬ 
tained from the magnetic Hamiltonian by integrating out the 
spatial dependence of the magnetization. The damping dyadic 
in the Thiele equation, Dj, can be obtained from the dissipa¬ 
tion function in the rigid core approximation, Vk(Xo), which is 
defined in the same way as in Eq. (5) but with the ansatz given 
in Eq. (18). For this system, it is more convenient to eval¬ 
uate the dyadic by performing the integration over all space 
after taking derivatives with respect to the core velocity. In 
other words, the dyadic can be obtained using the Lagrangian 
formulation by recognizing that 

£>t ■ X 0 = JJ- JJ dxdy JJ- [m, £>^ L (m) ihj ) . (21) 



FIG. 3. (Color online) In-plane magnetization fields for vortices and 
skyrmions. (a) Vector fields given by if>(x,y) in (18) for different 
values of q and c. (b) Vortex and (c) skyrmion for spin structure with 
c = 1 ,q = 1, where the arrows indicate the in-plane components 
{m xy ) and the color code gives the perpendicular component of the 
magnetization ( m z ). 


By using polar coordinates for the spatial coordinates, (x,y) = 
(r cosip, rsim/;), assuming translational invariance in the film 
plane, and integrating over ip, the damping dyadic is found to 
be 

Dj — -( ((XqDq + ct\D\ + ) I + CH2D2 

r \ 

( 22 ) 

where I is the 2x2 identity matrix and the dimensionless 
damping constants are defined as a\ = ?//r 2 , = >1 &r/ r c , and 

as = i] d 2 , in analogy with the domain wall case where the 
core radius r c plays the role here as the characteristic length 
scale. The coefficients D , depend on the core profile and are 
given by 


an 0 \ 

0 022 )’ 


D 0 =7T f 
Jo 

D\ = 2nr c 
D 2 — 2nr c 


*| r (Wc) J + Si, ' :e » 


rv 

Jo r 


r 

2 , • 2 . 


- (<9 r 0 o ) sin- 0 O 


(23) 

(24) 


0 r 

r°° 1 

I dr - (d r Oo) sin 0o (r (<3,-0o) cos 0o + sin 0 O ), 

Jo r 

(25) 


D 3 


-*r 


dr I r (<9,-0o) 2 cos 2 0 O 


sin“ 0o 


(26) 


where the expression for D 0 is a known result but D\, D 2 and 
Z >3 are new terms that arise from the nonlinear anisotropic 
damping due to RSOC. 

The coefficients a\\ and 022 are configuration-dependent 
and represent the chiral component of the Rashba-induced 
damping. For vortex-type spin textures (c - 1,3 and q = 1), 
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TABLE I. Coefficients a n and a 22 of the chiral damping term in 
Eq. (22) for different vortex/skyrmion charges q and chirality c. 


c 

0 

q = 1 

1 2 

3 

0 

q = 
1 

-1 

2 

3 

an 

1 

0 

-1 

0 

-1 

-1 

1 

1 

a 2 2 

1 

0 

-1 

0 

1 

-1 

-1 

1 


a 11 = 0-22 = 0 , which indicates that the a 2 term plays no role 
for such configurations. This is consistent with the result for 
Bloch domain walls discussed previously, since the vortex- 
type texture [Fig. 3(b)], particularly the vortex-type skyrmion 
[Fig. 3(c)], can be thought of as being analogous to a spin 
structure generated by a 2n revolution of a Bloch wall about 
an axis perpendicular to the film plane. The rigid core approx¬ 
imation implies that fluctuations about the ground state are ne¬ 
glected, which is akin to setting 54>{1) = 0 in Eq. ( 8 ). As such, 
no contribution from a 2 is expected for vortex-type textures. 
On the other hand, a finite contribution appears for hedgehog- 
type vortices and skyrmions (q = 1 ), where ci\\ = an = 1 
for c — 0 and an = a 22 = -1 for c — 2. This can be un¬ 
derstood with the same argument by noting that hedgehog- 
type textures can be generated by revolving Neel-type domain 
walls. A summary of these coefficients is given in Table I. 

For antivortices (q = -1), it is found that the coefficients 
a„ are nonzero for all winding numbers considered. We can 
understand this qualitatively by examining how the magneti¬ 
zation varies across the core along two orthogonal directions. 
For example, for c = 0, the variation along the x and y axes 
across the core are akin to two Neel-type walls of different 
chiralities, which results in nonvanishing contributions to an 
and a 2 2 but with opposite sign. The sign of these coefficients 
depends on how these axes are oriented in the film plane, as 
witnessed by the different chiralities c in Fig. 3. Such damping 
dynamics is therefore strongly anisotropic, which may have 
interesting consequences on the rotational motion of vortex- 
antivortex dipoles, for example, where the antivortex configu¬ 
ration oscillates between the different c values in time [41]. 

For vortex structures, we can provide numerical estimates 
of the different damping contributions ar,Z), by using the Usov 
ansatz for the vortex core magnetization. 


r > r c 

Let L represent the lateral system size. The coefficients Z), are 
then found to be Z> () = 7 r [2 + In (L/r c )], D\ = D 2 = 14zr/3, 
and /A = 71 [4/3 + In (L/r c )]. We note that for Do and /J 3 , 
the system size L and core radius r c appear as cutoffs for the 
divergent 1 /r term in the integral. By assuming parameters of 
oro = 0.1, q = 0.05 nm 2 , and a R = 0.1 eV nm, along with 
typical scales of r c = 10 nm and L - 1 /rtn, the damping 
terms can be evaluated numerically to be aoDo ~ 2.1, a\ I)\ * 
0.0073, a 2 D 2 ~ 0.19, and a 2 D 2 ~ 6.4. As for the domain 
walls, the Rashba term a^Dj, is the dominant contribution and 
is of the same order of magnitude as the linear damping term. 


cos 6q 


wink' the chiral term ow/T is an order of magnitude smaller 
and the nonlinear term a\D\ is negligible in comparison. 

For skyrmion configurations, a similar ansatz can be used 
for the core magnetization. 



r < r c 
r>r c 


(28) 


We note that this differs from the “linear” profiles discussed 
elsewhere [40], but the numerical differences are small and do 
not influence the qualitative features of the dynamics. The ad¬ 
vantage of the ansatz in Eq. (28) is that the integrals for D , 
have simple analytical expressions. Because spatial variations 
in the magnetization for a skyrmion are localized only to the 
core, in contrast to the circulating in-plane moments of vor¬ 
tices that extend across the entire system, the damping con¬ 
stants Dj have no explicit dependence on the system size. By 
using Eq. (28), we find D 0 = £>3 = 16zr/3, D\ = 4967 r/ 15 , and 
Z >2 = 52n/5. By using the same values of ao. //, and q-r as for 
the vortices in the preceding paragraph, we find aoDo ~ 1.7, 
a\D\ * 0.052, a 2 D 2 ~ 0.43, and 0 - 3^3 « 3.3. 


IV. DISCUSSION AND CONCLUDING REMARKS 

A clear consequence of the nonlinear anisotropic damp¬ 
ing introduced in Eq. (3) is that it provides a mechanism by 
which the overall damping constant, as extracted from domain 
wall experiments, for example, can differ from the value ob¬ 
tained using linear response methods such as ferromagnetic 
resonance [19]. However, the Rashba term can also affect the 
ferromagnetic linewidth in a nontrivial way. To see this, we 
consider the effect of the damping by evaluating the dissipa¬ 
tion function associated with a spin wave propagating in the 
plane of a perpendicularly magnetized system with an ampli¬ 
tude b(t) and wave vector ky. The spin wave can be expressed 
asm = [b(t) cos(k|| • ry ),b(t) sin(ky • ry), 1 ], which results in a 
dissipation function per unit volume of 

W sw = ^b(t) 2 (a 0 + a 3 + qb(tf ||ky|| 2 ), (29) 

where a term proportional to the chiral part qa^ spatially aver¬ 
ages out to zero. The Rashba contribution a 2 = qa\ leads to 
an overall increase in the damping for linear excitations and 
plays the same role as the usual Gilbert term ao in this ap¬ 
proximation, which allows us to assimilate the two terms as 
an effective FMR damping constant, ctfmr ~ a (> + a 3 - On 
the other hand, the nonlinear feedback term proportional to 
q is only important for large spin wave amplitudes and de¬ 
pends quadratically on the wave vector. This is consistent 
with recent experiments on permalloy films (in the absence of 
RSOC) in which the linear Gilbert damping was recovered in 
ferromagnetic resonance while nonlinear contributions were 
only seen for domain wall motion [19]. This result also sug¬ 
gests that the large damping constant in ultrathin Pt/Co/ALO, 
films as determined by similar time-resolved magneto-optical 
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microscopy experiments, where it is found that otfmr = 0 . 1 - 
0.3 [24], may partly be due to the RSOC mechanism described 
here (although dissipation resulting from spin pumping into 
the platinum underlayer is also likely to be important [42]). 
Incidentally, the nonlinear term rjb(t) 2 may provide a physi¬ 
cal basis for the phenomenological nonlinear damping model 
proposed in the context of spin-torque nano-oscillators [43], 

For vortices and skyrmions, the increase in the overall 
damping due to the Rashba term a 3 can have important con¬ 
sequences for their dynamics. The gyrotropic response to any 
force, as described by the Thiele equation in Eq. (19), depends 
on the overall strength of the damping term. This response 
can be characterized by a deflection angle, Oh , that describes 
the degree to which the resulting displacement is noncollinear 
with an applied force. This is analogous to a Hall effect. By 
neglecting the chiral term a^A, the deflection or Hall angle 
can be deduced from Eq. (19) to be 


a 0 D 0 + a\D\ + a^A ' 

where Go = 2n for vortices and Go = 47 r for skyrmions. Con¬ 
sider the skyrmion profile and the magnetic parameters dis¬ 
cussed in Section III. With only the linear Gilbert damping 
term (oroA)) the Hall angle is found to be Oh — 82.3°, which 
underlies the largely gyrotropic nature of the dynamics. If 
the full nonlinear damping is taken into account [Eq. (30)], 
we find Oh = 68.3°, which represents a significant reduction 
in the Hall effect and a greater Newtonian response to an ap¬ 
plied force. Aside from a quantitative increase in the overall 


damping, the presence of the nonlinear terms can therefore af¬ 
fect the dynamics qualitatively. Such considerations may be 
important for interpreting current-driven skyrmion dynamics 
in racetrack geometries, where the interplay between edge re¬ 
pulsion and spin torques is crucial for determining skyrmion 
trajectories [44, 45]. 

Finally, we conclude by commenting on the relevance of 
the chiral-dependent component of the damping term, a 2 . It 
has been shown theoretically that the Rashba spin-orbit cou¬ 
pling leading to Eq. (3) also gives rise to an effective chiral 
interaction of the Dzyaloshinskii-Moriya form [34]. This in¬ 
teraction is equivalent to the interface-driven form considered 
earlier, which favors monochiral Neel wall structures in ul- 
trathin films with perpendicular magnetic anisotropy. Within 
this picture, a sufficiently strong Rashba interaction should 
only favor domain wall or skyrmion spin textures with one 
given chirality as determined by the induced Dzyaloshinskii- 
Moriya interaction. So while some non-negligible differences 
in the chiral damping between vortices and skyrmions of dif¬ 
ferent chiralities were found, probing the dynamics of solitons 
with distinct chiralities may be very difficult to achieve exper¬ 
imentally in material systems of interest. 

ACKNOWLEDGMENTS 

The author acknowledges fruitful discussions with P. Borys, 
J.-Y. Chauleau, and F. Garcia-Sanchez. This work was par¬ 
tially supported by the Agence Nationale de la Recherche un¬ 
der Contract No. ANR-ll-BS 10-003 (NanoSWITI) and No. 
ANR-14-CE26-0012 (ULTRASKY). 


[1] L. Landau and E. Lifshitz, Phys. Z. Sowjet 8, 153 (1935). 

[2] M. Sparks, Ferromagnetic-relaxation theory (McGraw-Hill, 
New York, 1964). 

[3] T. Gilbert, IEEE. Trans. Magn. 40, 3443 (2004). 

[4] Y. Tserkovnyak, A. Brataas, G. E. W. Bauer, and B. Halperin, 
Rev. Mod. Phys. 77, 1375 (2005). 

[5] Y. Tserkovnyak, A. Brataas, and G. E. W. Bauer, Phys. Rev. 
Lett. 88, 117601 (2002). 

[6] D. Mills, Phys. Rev. B 68, 014419 (2003). 

[7] H. Hurdequint, J. Magn. Magn. Mater. 93, 336 (1991). 

[8] S. Mizukami, Y. Ando, and T. Miyazaki, Jpn. J. Appl. Phys. 
40, 580 (2001). 

[9] S. Mizukami, Y. Ando, and T. Miyazaki, Phys. Rev. B 66, 
104413 (2002). 

[10] H. Hurdequint and M. Malouche, J. Magn. Magn. Mater. 93, 
276(1991). 

[11] R. Urban, G. Woltersdorf, and B. Heinrich, Phys. Rev. Lett. 87, 
217204 (2001). 

[12] J.-V. Kim and C. Chappert, J. Magn. Magn. Mater. 286, 56 
(2005). 

[13] X. Joyeux, T. Devolder, J.-V. Kim, Y. Gomez De La Torre, 
S. Eimer, and C. Chappert, J. Appl. Phys. 110, 063915 (2011). 

[14] S. Zhang and S. S.-L. Zhang, Phys. Rev. Lett. 102, 086601 
(2009). 

[15] L. Berger, J. Appl. Phys. 55, 1954 (1984). 


[16] S. Zhang and Z. Li, Phys. Rev. Lett. 93, 127204 (2004). 

[17] D. Claudio-Gonzalez, A. Thiaville, and J. Miltat, Phys. Rev. 
Lett. 108, 227208 (2012). 

[18] A. Manchon, W. S. Kim, and K.-J. Lee, “Role of Spin Diffusion 
in Current-Induced Domain Wall Motion,” arXiv: 1110.3487 
[cond-mat:mes-hall]. 

[19] T. Weindler, H. G. Bauer, R. Islinger, B. Boehm, J. Y. Chauleau, 
and C. H. Back, Phys. Rev. Lett. 113, 237204 (2014). 

[20] K.-W. Kim, J.-H. Moon, K.-J. Lee, and H.-W. Lee, Phys. Rev. 
Lett. 108,217202(2012). 

[21] H.-B. Braun, Adv. Phys. 61, 1 (2012). 

[22] C. Burrowes, N. Vernier, J.-P. Adam, L. Herrera-Diez, K. Gar¬ 
cia, I. Barisic, G. Agnus, S. Eimer, J.-V. Kim, T. Devolder, 
A. Lamperti, R. Mantovan, B. Ockert, E. E. Fullerton, and 
D. Ravelosona, Appl. Phys. Lett. 103, 182401 (2013). 

[23] T. Devolder, P. H. Ducrot, J.-P. Adam, I. Barisic, N. Vernier, J.- 
V. Kim, B. Ockert, and D. Ravelosona, Appl. Phys. Lett. 102, 
022407 (2013). 

[24] A. J. Schellekens, L. Deen, D. Wang, J. T. Kohlhepp, H. J. M. 
Swagten, and B. Koopmans, Appl. Phys. Lett. 102, 082405 
(2013). 

[25] A. Thiaville, J. Garcia, and J. Miltat, J. Magn. Magn. Mater. 
242, 1061 (2002). 

[26] Y. Le Maho, J.-V. Kim, and G. Tatara, Phys. Rev. B 79, 174404 
(2009). 




[27] A. N. Bogdanov and U. K. RoBler, Phys. Rev. Lett. 87, 037203 
( 2001 ). 

[28] A. Thiaville, S. Rohart, E. Jue, V. Cros, and A. Fert, Europhys. 
Lett. 100. 57002 (2012). 

[29] A. Fert and P. M. Levy, Phys. Rev. Lett. 44, 1538 (1980). 

[30] A. Fert, Mat. Sci. Forum 59-60, 439 (1990). 

[31] M. Bode, M. Heide, K. Von Bergmann, P. Ferriani, S. Heinze, 
G. Bihlmayer, A. Kubetzka, O. Pietzsch, S. Bliigel, and 

R. Wiesendanger, Nature 447, 190 (2007). 

[32] M. Heide, G. Bihlmayer, and S. Bliigel, Phys. Rev. B 78, 
140403 (2008). 

[33] H. Imamura, P. Bruno, and Y. Utsurni, Phys. Rev. B 69, 121303 
(2004). 

[34] K.-W. Kim, H.-W. Lee, K.-J. Lee, and M. D. Stiles, Phys. Rev. 
Lett. Ill, 216601 (2013). 

[35] M. Belmeguenai, J.-P. Adam, Y. Roussigne, S. Eimer, T. De- 
volder, J.-V. Kim, S. M. Cherif, A. Stashkevich, and A. Thiav¬ 
ille, Phys. Rev. B 91, 180405(R) (2015). 

[36] I. M. Miron, T. Moore, H. Szambolics, L. D. Buda-Prejbeanu, 

S. Auffret, B. Rodmacq, S. Pizzini, J. Vogel, M. Bonfim, 


A. Schuhl, and G. Gaudin, Nat. Mater. 10, 419 (2011). 

[37] J. P. Tetienne, T. Hingant, L. J. Martinez, S. Rohart, A. Thiav¬ 
ille, L. H. Diez, K. Garcia, J.-P. Adam, J.-V. Kim. J. F. Roch, 

I. M. Miron, G. Gaudin, L. Vila, B. Ocker, D. Ravelosona, and 
V. Jacques, Nat. Commun. 6, 6733 (2015). 

[38] E. Feldtkeller and H. Thomas, Z. Phys. B 4, 8 (1965). 

[39] Y. Gaididei, V. P. Kravchuk, and D. D. Sheka, Int. J. Quant. 
Chem. 110, 83 (2010). 

[40] N. S. Kiselev, A. N. Bogdanov, R. Schafer, and U. K. RoBler, 

J. Phys. D: Appl. Phys. 44, 392001 (2011). 

[41] G. Finocchio, O. Ozatay, L. Torres. R. A. Buhrman, D. C. 
Ralph, and B. Azzerboni, Phys. Rev. B 78, 174408 (2008). 

[42] J.-M. L. Beaujour, J. H. Lee, A. D. Kent, K. Krycka, and C.-C. 
Kao, Phys. Rev. B 74, 214405 (2006). 

[43] V. Tiberkevich and A. Slavin, Phys. Rev. B 75, 014440 (2007). 

[44] A. Fert, V. Cros, and J. Sampaio, Nat. Nanotech. 8, 152 (2013). 

[45] J. Sampaio, V. Cros, S. Rohart, A. Thiaville, and A. Fert, Nat. 
Nanotech. 8, 839 (2013). 



